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The nonlinear equations that represent the behavior of valve-controlled 
hydraulic servomechanisms are derived, and the assumptions necessary for 
their linearization are discussed. Solutions of the nonlinear equations ob- 
tained by analog computation arc compared with solutions of the linear equa- 
tions. Attention is directed to the influence of the hydraulic parameters on 
the nonlinear closed-loop system behavior. 

I. INTRODUCTION 

Since the development of hydraulic control valves such as that em- 
ployed in the Nike missile,' emphasis has been given to the analysis 
of hydraulic phenomena in valves! systems,--^* with much of the litera- 
ture having been devoted to hydrauhc component design. In this paper, 
the nonlinear closed-loop performance is given major emphasis; the effect 
of hydraulic parametei- variations on the closed-loop frefjuency and 
transient responses is examined by linear and nonlinear methods. 

The basic servomechanism under consideration in this study, as shown 
in Fig. 1, consists of a summing device, an amplifier, a flow source and 
control system, a hydraulic actuator (or motor) and a load. In the 
following sections, the nonlinear differential equations which represent 
the l>ehavior of this closed-loop system are derived, a linear and an 
incremental-linear representation are discussed and solutions of the 
nonlinear equations obtained by analog computation are compared with 
linear solutions. 

II. MOTOR AND I.OAn ANALYSIS 

In this section, emphasis is on the derivation and validity of the 
equations used to represent the l)ehavior of the actuator and load; the 
mechanization of the flow source and its method of control will be 
discussed hi detail in the following section. 
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Fig. 1 — Basic hydraulic positional servomechaniam. 

In Fig. 2, a pi.ston-type actuator is shown connectoci hydraulically 
to a flow source and mechanically to a load. The initial step in the 
analysis is to relate the flows Qi and Qa to the dependent variable x, 
the piston displacement from the center position. 

The instantaneous volumes between the piston and two arl>itrary 
sections in the lines leading to the cylinder are designated Vji and 
Vt2 (in cubic inches), the numerical subscripts indicating a particular 
side of the piston. Qi and Q^ (in cubic inches per second) represent the 
flows from the source and depend upon the som'ce mechanization; they 
are considered positive in the directions indicated in the diagram. Ql , 
the leakage flow past the piston, is also shown in the assumed positive 
direction. 

A control volume is chosen so that it coincndew with the volume Vn 
(where Vn is a function of time) and the cfiuation for the conservation 
of mass flow is written for this volume. This relationship states that the 
rate of mass accumulation in the (control volume is e(iual to the net 
rate of mass flow into the volume. The net rate of mass flow into Vn is 
given by 

net rate of mass flow into Vn = p(Qi — Qi.), (1) 

where p is the mass density (in pound-seconds^ per inch^) of the fluid 
and Qi. is the flow out of Vn ■ 
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Fi}?. 2 — Variable diw placement actuator. 
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In writing (1), it is assumed that the mass density, p, is uniform 
throughout the control vohmie; i.e., it is assumed that p is a function of 
time only, and uot a function of position in the volume. Since p is de- 
pendent upon the instantaneous pre.ssure, pi , existing in Fn , this 
assumes that pi is miiform throughout Vn ■ The justification for this 
assumption is based on the calculation of the velocity of propagation 
of a longitudinal compression wave in the fluid. In general, the velocity 
of propagation, for adiabatic conditions, is given by 

v.= ^l, (2) 

where /3 is the adiabatic luilk modulus of compression of the fluid. For 
a typical oil, Vp is approximately 50,000 inches per second. If the largest 
linear dimension of the volume T'ti is small, the pressure wave will 
complete many cycles in a short time, and the perturbation will be 
rapidly attenuated. For valve-controlled high-performance systems, 
the volumes are small; imder the assumption that the largest dimension 
is one inch, the time of travel is 0.02 millisecond. Since this value is 
small compared to system time constants of the order of five milliseconds 
or more, the pressure can be assumed uniform throughout the volume. 
A calculation of the frequency of oscillation at which nonuniform pres- 
sure distribution becomes important shows that it is much higher than 
the frequencies of interest: 10,000 cps versus 200 cps. 

The rate of mass accumulation in the volume Ft-i is given by 

rate of mass accumulation in Vri = .. = pFj-i + Fnp. (3) 

at 

Equating (1) to (3) and solving for Qi yields 

Qi - Ql + Fti + ^ p. (4) 

p 

The adiabatic bulk modulus of compression, /3, of the fluid is defined as 

^-^^' (5) 



(?)^ 



where /3 is assumed constant. The ehmination of p from (4) by use of 
(5) results in 
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In the same manner, application of the equation for the conserva- 
tion of mass flow to the volume Fra yields 

Qi = -Ql+ Vt^+^V2. (7) 

Experimental tests show that the leakage flow is proportional to the 
pressure differential across the motor (laminar flow), so that 

Qj. = L,„{p, - p2), (8) 

where Lm is the leakage coefficient (in inches per pound-second) of 
the actuator. In addition, define 

y, = Yil^J^^ (9) 

where Vt is a constant, so that 

Vr, = Vt-\- Ax, 

(10) 
Vt2 = Vt — Ax, 

where A is the cross-sectional area (in square inches) of the piston and 
X (in inches) is the piston displacement measured from the center 
position. Equations (6) and (7) are thus 

Qi = L^(V. - P2) -\-Ax-\- ^^"^^"^^ pi (11) 

and 

Q2 = -Ljp, - P2) ~Ax-\- ^^" ~ ^""^ P2 . (12) 

These flow equations have been developed for the linear piston-type 
actuator, but the same equations are valid for vane motors. 

For a fixed-stroke axial-piston rotary motor, the control volume Vn 
is a discontinuous function of time, since, as the cylinder block rotates, 
the individual cylinders transfer from one side of the motor to the 
other. Since the volume of one cyfinder is small compared to the total 
volume on one side of the motor, the volume variation due to this 
discontinuity may be neglected without serious error. Thus, the control 
volume V-n is essentially constant, so that 

Fn - Vt2= Vt, (13) 

where Vt is constant and is equal to one-half of the volume in the 
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system. With this assumption, the flow equations for the fixed-stroke 
axial-piston motor can be written 

Qi = L^Pi - P2) + VJ, + y Pi (14) 

and 

Q2 - -L,„(pi-p2)-FJo+ yP2, (15) 

where ^o is the motor shaft angular rate and V,„ (in cubic inches per 
radian) is the fluid displacement per unit rotation of the motor shaft. 
Since Vt and do not occur separately in (14) and (15), it will be 
convenient to define a "compliance coefficient" Kc as 

/c = ^^ (16) 

where Ka has the units of inches" per pound. 

The remaining discussion will be c(tnccrned with a rotational system, 
for which (14) and (15) have been developed; the same relationships 
will be vafid for a translational system having small displacement x, 
with appropriate changes in the definitions of the parameters. 

In addition to the flow equations, two torque equations can be derived. 
The first is an energy relationship that equates the w^ork done by the 
forces on the motor during a rotational of ^o radians to the work output 
from the motor shaft. The work input is (pi — p^) V.A and is the flow 
work commonly encountered in the Bernoulh equation. The corre- 
sponding work output is Td^ , where T is the opposing toniue. Since 
these two expressions for work must be equal, there results 

T = (pi- P2) V,.. (17) 

Another torque equation is obtained from Newton's Second Law^ of 
Motion. In general, the load may consist of inertia, damping and fric- 
tion toniues, and disturbing torques. Thus, the following equation may 

be written: 

T = J9,+ To, (18) 

where ./ (in pound-inch-.seconds") is the total inertia (including that 
of the motor and fluid) referred to the motor shaft, and Tq (in pound- 
inches) is the total friction and disturbing torque acting on the motor 
shaft. 
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The elimination of T from (17) and (18) results in 

(pi - p.) V,. - Jk+ n. (19) 

A restriction must be placed on the values allowable for the pressures 
Pi and p2 ■ If absolute pressure units are employed, these pressures 
must always be equal to or greater than zero. A more accurate repre- 
sentation would be obtained if the vaporization pressure were con- 
sidered as the limithig value, but, in view of the fact that the differential 
pressures in the system are normally very large, this degree of refine- 
ment is not warranted. A good approximation is 

pi g and Pa g 0, (20) 

where it is understood that absolute pressure units are employed. 

In addition to the equations derived, the usual expression for the 
position error e in terms of the input angle dt and the output angle 6o 
for a Hcrvomechanism having unity feedback (as in Fig. 1) is given by 

e^ di- e^. (21) 

The equations which have been derived in this section and which 
apply to the axial-piston rotary motor are summarized below: 

Qi = L,„(pi - p.) + F Jo + K4, , 

Q, - -L,„(pi - P2) ~ v., A + KcP2, 

(pi - P2) V,. = ,19, + To , (22) 

Pi ^ and p2 ^ 0, 



e = Ui ~ 



The units employed in these equations are given in Table I. 

The expressions for Qi and Q2 , the flows from the controlled source, 

are discussed in the next section. 

III. CONTROL-VALVE ANALYSIS 

In Section II, the equations relating the flows (into the motor) to 
the dynamic state of the system were derived; the expressions for the 
flows Qi and Q2 were not specified. In this section, these quantities are 
discussed for the particular case of a valve-controlled system and 
analytical expressions relating flow to error signal and pressure are 
obtained. 

The schematic diagram of a typical valve configuration is given in 
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Table I — Definitions of Symbols; Units 



Symbol 


Definition 


Units 


Q 


flow 


in^/sec 


L^ 


motor leakage coefficient 


in'' /lb-sec 


V 


pressure 


Ib/in^ 


F™ 


motor displacement 


in^/radian 


di 


position angle input 


radians 


60 


position angle output 


radians 


e 


position angular error 


radians 


Vt 


one-half total trapped volume 


in'' 


P 


bulk modulus of compression 


lb/in= 


J 


total inertia 


Ib-iu-sec'^ 


Kc 


compliance coefficient 


inVlb 


To 


viscous, friction and disturbing tovqiie 


Ib-in 


p 


fluid mass density 


Ib-see^/in* 



Fig. 3. The type of actuator is not important in this discussion and, for 
simplicity, it is pictured as a translational piston. 

A source having a pressure p^ supplies fluid to the valve as shown 
and the mam spool controls the direction and magnitude of this flow 
to the motor. Fluid is returned to a sump at pressure pd . As in Section 
II, the pressure on each side of the actuator piston is designated as p 
with the appropriate numerical subscript. 

The main spool controls the flow by means of four orifices Oi , O-i , 
El and E^ , and the spool position, in turn, is controlled by a transducer. 
The configuration of the transducer varies considerably in current 
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Fig. 3 — Orifice flow conventions. 
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production models and may have the form of either a torque motor or 
a hydraulic preamplifier preceded by a torque motor or solenoid. If a 
hydraulic preamplifier is employed, the entire valve is designated as a 
"two-stage" valve. In the analysis an ideal spool-positioning mechanism 
will be assumed; i.e., there will be only one position of the valve spool 
corresponding to a given error signal, e. 

The following nodal equations are obtained from Fig. 3: 

Qi — Q()2 — Qez ■ 

Here, as in Section IT, the flows are chosen positive in the directions 
indicated in the figure. 

It is necessary to express the flows through the orifices in terms of 
the pressures and the orifice openings. Consider first the general orifice 
equation; application of Bernoulli's equation to the case of flow through 
an orifice of area A results in the relation 



Q = 



AC,C\ , /2Ap 



i/'-.'(0 



P 



(24) 



where : 



Q = flow, 

A = orifice area, 

Cc = orifice contraction coefficient, 

C„ — velocity coefficient, 

Ai = upstream line area, 

Ap = pressure differential across the orifice, 

p = mass density of the fluid. 

For appUcation to a valve orifice, note that A/Ai is small compared to 
unity and that Cc is less than unity. It follows, therefore, that a good 
approximation is obtained by 

Q = AC.C.^^. (25) 

The usual procedure in hydrauhcs is to define a discharge coefficient, 
C, as 

C = CcC, , (26) 

so that (25) becomes 

Q^AC y^^. (27) 
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This equation has been derived for the steady-state flow through an 
orifice; m the analysis, it will be assumed that the same relationship is 
valid for the dynamic state. It is not necessary to assume a value for 
the discharge coefficient C, since it will be included in an over-all gain 
constant. 

In the present application, the orifice area, A, is proportional to the 
displacement of the valve spool (assuming rectangular orifices) and 
the displacement is proportional to the electrical activating signal, Ke, 
received by the positioning mechanism (see Fig. 3). It follows that A 
is proportional to Kf, and (27) may be written as 

Q = (Ke)C,,y^^, (28) 

where Co is a. new coefficient that is proportional to the area of the 
orifice per radian of inp\it signal. If the mass density, p, of the fluid is 
assumed constant in this expression, the following equation may be 
written : 

Q - (K,e) VAp, (29) 

where 

7^, = ^^^' (30) 

is a constant for a given system. Kb has the units of inches* per pound^'^- 
second-radians. 

The relationship of (29) is indicated graphically in Fig. 4. In this 
figure, the equation has been normalized with respect to the three 
maximum values Qmux , {Kbe),na^ and Ap,,,,,.^ . Since there is a limit to 
the magnitude of the spool displacement, a spool-displacement satura- 
tion region is indicated. (Spool-displacement saturation occurs when 
the electrical error signal becomes greater than that corresponding to 
the maximum valve spool-displacement; i.e., the error signal demands 
a spool position that is physically impossible.) 

Equation (29) implies that the flow through the orifice is zero for 
zero error e. This is not generally true for the orifices in most valves, 
and in fact, this condition would be very difficult to obtain. Equation 
(29) must be specialized for each orifice. 

There arc three general types of valves, which can be classified accord- 
ing to the flow conditions at zero signal input e; these are: (a) the open- 
center valve, (b) the critical-center valve and (c) the closed-center 
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Fig. 4 — Orifice flow characteristics. 
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Fig. 5^ (a) Neittriil position for tlic three basic valve types; (b) typical 
characteristics for the three basic valve types. 
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Fig. G — Flow cliaraeteristics for the orifices in au open-ceuter valve. 

valve. Fig. 5(a) indicates the zero error position of the spool for each 
of these types for the orifice Oi . In the open-center valve, flow passes 
through the orifice for the condition e = 0, i.e., for the spool in the 
neutral position. The critical-center valve allows no flow to pass in the 
neutral position; however, a slight positive displacement opens the 
orifice 0\ . The closed-center valve has a dead zone in that a relatively 
lart^e tlisplacement is required to open the orifice from the neutral 
po.sition. 

Kig. 5(b) shows a typical differential i.sobar (corresponding to those 
of l'"ig. 4) for each type of vahe. The origin is translated according to 
the neutral spool position. 

It is now possible to express the individual orifice flows as shown in 
Fig. 3 with the aid of Fig. 6 and (29). Fig. 6 shows a typical differential 
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isobar for each of the four orifices Oi , O2 , -Ei and & , the case shown 
representhig the open-center type valve. (The equations to be derived 
will be applicable to all three valve types, depending upon the choice 
of the neutral poHition error-flow constants.) It should be noted that 

the direction of flow is dependent on the sign of Ap for each orifice. In 
general, the equations for the flows arc as given by (31). In these equa- 
tions €01 , €02 , ^Ei , e£2 fire positive for an open-center valve, en,„x is 
taken as positive and sgn denotes the signum (sign) function: 

'Kb\(:,n^.K + em I V\Ps — pi I Hgn (p, - pi) e > e„mx 
Qoi ^ Kb\e -\- €o\ \-\/\ps — Pi\ Sgn (p, — pi) -601 < e < e^^^^ 

_0 € < — €01, 

e > eo2 

Q02 ^ Kh\€ - eo2 ] Vl ps - P2 I sgn (p, — P2) —e„uix <€ < eo2 

_Kb I emax — «02 | V | Ps " ^2 | SgU (p« — ^2) e < — Cmiix , 

_ (31 J 

e > e^i 

Qfi-i = Kb\€ — cei j VI pi - pd] sgn (pi - Pd) -e,nax < € < CeI 

_Kb I Cmnx — tiJl I Vl Pl — P<Z I Sgn (pi — Pd) f < — Cmax , 

Xh [ tmiiK + €E2 I V|P2 — Pd| Sgn (p2 — Pj) e > Cmax 

Qfi3 ^ -Kft I « + «B2 I V| P2 - Pd I Sgn (p2 - Pd) —fE-i < € < Cmax 

e < — (£2 ; 

Simplification of (31) is obtained by assuming that the sump pressure 

Table II — Definitions of Symbols; Units 



Symbol 


Definition 


Units 


Q 


flow 


inVsec 


A 


orifice area 


in^ 


Ai 


upstream line area 


in^ 


c. 


orifice contraction coefficient 


— ■ 


c. 


orifice velocity coefficient 


■ — 


C 


orifice discharge coefficient 


— 


Ct 


area per radian input 


in^/rad 


t 


angular actuating signal 


radians 


K 


dimension! ess gain constant 


— 


Kb 


over-all gain constant 


inVlb'^'i-sec-rad 


P 


pressure 


Ib/in^ 


Ap 


differential pressure 


Ib/in^ 


P 


fluid mass density 


Ib-sec^/in^ 



VALVE-CONTEOLLED HYDR^VULIC SERVOMECHANISMS 



1525 



p,{ is zero psi absolute rather than atmospheric pressure; this is a good 
appi'oximatioii since the pressure differentials are usually several hundred 
psi. This assumption eliminates consideration of the shaded regions in 
Fig. for the orifices ^i and E-i . Further simphfication is obtained if 
the valve is assumed to be perfectly symmetrical; i.e., 

and if it is assumed that the maximum spool po.sition is never attained. 
With these assumptions, (31) simplify to (32): 

/U I c + eo I V\ P» - Vi\ sgn {ps -pi) €> -fo 

e < -eo, 

e > eo 

_Kh\c - eolVlps - P2I sgn {p, - p.) e < to, 

"0 e> €0 

_Kb I e — to I Vpi e < CO , 

Kh\t + f u I ■\/v'i e > — eo 

e < -Co. 



Q02 = 
Qei = 
Qei = 



(32) 



These equations, together with (23) and (22), complete the prelimi- 
nary analysis for the valve-controlled servo mechanism. Definitions of 
the symbols, together with a consistent set of units, are given in Table II. 

IV. LINEAR AND INCREMENTAL-LIXEAH ANALYSIS 

The linearization of the equations representing the motor as given 
by (22) results in (for a pure iuertial load): 

(h = LJp, - p.) + F,A + /Cpi , 
(h = -LJpi - P2) - r,A + /wp2, 

(pi — ps) Fm = J6n, 



(33) 



f = 0, — 



The restrictions on the \'alnes of pi and P'> and their derivatives are 
not applicable to a linear theory and have been omitted; in addition, 
the friction and output disturbing torques have been assumed to be zero. 
For the case of a symmetrical valve having characteristics as given 
by (32), the flows through the orifices are functions of the gain constant, 
Kb , the error signal, e, the open-center constant, fo , and the respective 
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orifice prewsure drops. If the region of applicability is restricted to those 
cases in which -p, > pi and p^ > pi {pressures not limited) and in which 

I e I < I CinHx I, (32) maj' be written: 

Qoi = /u(e + eo)V'p^ — Pi 1*^1' * > ~*o ' 

Qo2 ^ — A"b(e — eo)V'p« — P2 ft)!' e < «u, 

Qm ^ — /u(« — *(.)\/pi foi' e < f , 

Qk2 ^ Kb(e + fo)A/p^ for « > — «o • 

From (23) and (34) it can be seen that there are three distinct ranges 
of € that must be considered, the magnitude of the range depending 
npon the value of eo : 
Region A (| e | ^ en) : 

Qi ^ Kb(€ + eo^Vp. — Pi + A'b^^ — eo)-\/pi ) 
C^s - -Kb(e - eo)\/p« - P-J - ^^''■(^ + eo)-s/pa ■ 
Region B-\- (e ^ cu) : 



(35) 



f^i - /u(e + eo)\/ps - Pi ) 
^^2 = -/u(e + 6u)a/p2 . 
Region B— (e ^ to) : 

Qi = Kb(€ - eo)-\/pi, 

Qi ^ —Kbi^ — eo)\/p« — P2 



(36) 



(37) 



Expanding \/ps — p f^^cl \/p about the steady-state value Pi/2 
yields 
Region A : 



(38) 



ai-C2-2/^.(2.^|"-.„^;7^^). 

Region B+ : 

0, - 0, = lu U + ..) (2 |/| - ^^') . (39) 

Region B— : 

Qi-Q. = lUe - eo) (2 j/'l + ^^7^')> (40) 

where the higher-order terms in [(pi,/2) — p] have been neglected. This 
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further restricts llie region of ^-alidity of our analysis to those cases 
where the pressure differential across the load is not large. 
For the region A, combination of (33) and (38) results in 

4;u y'l e = ^:^ h + f^^ (l, + ^) h + 2Y.A , (41 ) 
and the open-loop transfer function becomes 



gu(.s-) KJ y 2 

:(S) 



G,{s) ^''-^ =— ^ /'^ l_\ ^77^, (42) 






or 



GAs) = -^-r = , ^ ^y ^, 7-^ =,, 143, 



t^n/ L\"'i/ " Vt*)-./ J 



where 



^.- K,/</#7, (44; 



FtJ 



1 



:^-+^T]i/t' <^^' 



2A', /p. 



(46) 



Here, u„ is the undamped natural frequency of the system, f„ is the 
clamping ratio and oj,, is the velocity gain constant. (The suliscript "a" 
indicates the region A.) It is interesting to note that t<i , the damping 
ratio, is the sum of the motor damping (motor leakage) and a term 
related to the steady-state flow through the valve. The term, Kbea/-\/'2p, , 
contributes the major damping to the system. In the quiescent state, 
siuee 6 = and pi = p-i ^ p.. 2, it follows from (34) that 



(Qoi), = (Q02). - (Qei)s = {Qe-2). = K,eo /j/^, 



(47) 



where the subscript "s" indicates quiescent values. Xow designate 

a = {Qoi).+ (Q02)., (48) 

where Q, is the quiescent (i.e. e =^ and 0^0) total flow from the 
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source, so that 



Q, = 2if.6o/j/|- (49) 

An effective leakage coefficient is given by 

and therefore (45) becomes 

Now consider the region B+, where e ^ eo , The combination of 

(33) and (39) results in 

2/C.y^l (e 4- .o) -^/o +^[Z-. +^^75^]^"+ 2^"'^«- <52) 

An approximate "incremental linear" transfer function may be obtained 
for this region by making the substitutions 

e - e* + At, 

(53) 

du = 6,* + A^o , 

where both the starred and tlie incremental symbols are considered as 
functions of time. The incremental variables are assumed small, so that 
their products may be neglected. In this manner, a linear equation in 
the incremental quantities can be obtained if the equation defining the 
starred variables [(52) with the symbols .starred] is subtracted from 
that obtained by the substitutions indicated previously. Then, if the 
resultant equation, which is hnear in the incremental variables, con- 
tains any starred quantities, these can be considered to be varying 
slowly with time — that is, essentially constant when compared with 
the incremental variations with time. Therefore, a quasilincar incre- 
mental transfer function can be obtained. In the present case, the 
incremental transfer function valid for a small constant-acceleration, 
do , will be derived. This will be used to obtain the incremental transfer 
function for constant-velocity operation. 

Thus, the analysis is initiated by obtaining the linear equation in 
the incremental quantities as previously outlined; this equation is 
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Noting that a similar equation may be obtained for the region B — , 
and denoting both regions l)y B, the following transfer function is 
obtained for operation alxmt a constant acceleration ao : 

^ . , Ado(s) o),, (r:^) 

"'" (£)[(£) -4£)-0 



wliere 



- V„ 



Z' ''"' 



-":i/i(^-fel)- <««> 

For an inertia load (no viscous or coulomb friction), (58) can be written: 

where Ap* is the constant-differential pressure acting on the motor. 

The preceding equations give an approximate solution for the case 
in which the acceleration ao is small; i.e., «* varies slowly with time. 
For this case, the incremental damping ratio ft is appreciably increased 
over that given by (45) ft)r the region A. The incremental gain constant 
016 is less than the gain constant in the A region as defined in (40); for 
Ap* equal to p, , the gain is down 12 db from that of the A region. 

For operation about a constant velocity on , w„ and fb remain as in 
(56) and (57), but the gain constant becomes 

Kb Jp. ^gO) 



' 111 F — 

In this case, the gain is down db from that gi\-cn for the A region. 
For constant-speed operation, ^o* = w,i, the error is given by 

so that the effective leakage coefficient becomes: 

U = L.., + ^^^''''^*j + ^"^ - L. + ^^ (62) 

2\/2p. ^P^ 
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Tlnis, since the "di.spl.acemeiit" motor flow is given by 

Q.„ - V,A = V«Mi , (63) 

the effective leakage coefficient is 

L, = L,„ + |^. (64) 

In many cases, L„, is small compared to Q™ so that the incremental 
damping is primarily a function of the motor speed. 

The velocity-lag error, e, (steady-state position error under the 
conditions ^o = o-h where w, is a constant rate input), is as follows for 
the two regions: 

ffegion A : 



Region B : 






eo, (06) 



where it is understood that | e^, \ is greater than co . The velocity lag 
error for the B region is thus approximately twice that predicted by 
(05) if Co is small compared to the lag error. 

Some qualitative information on the nature of the system performance 
can be obtained by comparison of (43) and (55). In the region A (i.e., 
j E { < Co), the system is essentially hnear for small pressure differentials. 
In the region B (I e | > tu), the system is nonlinear even though the 
pressiu'e differentials are assumed small. In this region, the behavior 
is "amplitude sensitive"; as the error amplitude increases, the incre- 
mental gain decreases. The incremental damping increases with increas- 
ing error in the B region and is considerably greater than the damping 
in the A region; the incremental damping is proportional to the total 
flow through the \^alve. The \'elocity-lag error is, of course, greater in 
the B region. 

The transfer function Ga{s) as given by (43) was derived for a fully 
symmetrical valve; the analysis of an unsymmetrieal valve shows that 
the liasic form of the transfer function is similar to that for the sym- 
metrical valve. For the unsymmetrieal valve, however, the damping 
ratio and gain expressions differ from those of (45) and (40). The damp- 
ing ratio is 



VALVE-CONTROLLED HYDRAULIC SERVOMECHANISMS 1531 






Avhere p, is the {luicscent value of the pressures pi and p2 . For p, = 
p„/2, this expression is identical to that given Ijy (45) ; for p., greater or 

less than ps/2, thcdiimpingis greater than that given by (45) (assuming 
the same quiescent flow Q^). The gain constant is: 

Oia = J^ (Vpi + Vp- - Pa)y (68) 

and comparison with (4fi) shows that the two give the same solution 
for p,, = pa/'2. For pq greater oi- less than p,/2, the gain of the unsym- 
metricnl valve is less than that of the symmetrical valve. 

V. SOLUTIONS OF THE LINEARIZED EQUATIONS 

The equations representing the valve-controlled servomechanism 
were linearized in the last section, and it was found that the open- 
loop transfer function had the following general form: 

G(s) = ^ ^ , . ^. ,. "" —^ =i, (69) 

e(s) 



(£)[(£)• -t.)-]^ 



where w,, is the velocity gain constant, aj„ the natural resonant frequency 
and i" the dimensionless damping ratio. If this is considered as a fre- 
quency function, the resulting open-loop attenuation and phase vary 
as shown in Figs. 7 and 8. In these figures, the damping ratio, f, has 
been taken as a parameter. 

Equation ((iO), when solved for the closed-loop function, results in 



ft,(.0 ^ 

his) ^-^^ 



0)0 






\w„/ \U},J 



(70) 



The relationships for the closed-loop opei'ation are exhibited graphi- 
cally in Figs. 1) through 12. In Fig. 9, the gain margin is shown as a 
function of the peak attenuation, i]fp (the maximum value of | V^< |), 
for the range of values of interest. In general, the gain margin decreases 
with increasing peak magnitucU^ and is less for the lightly damped cases. 
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Fig. 7 — Open-loop attenuation of G{,s). 

The phase margin variation with the peak magnitude is shown in Fig. 
10; the smaller values of f have the largest phase margins. 

Fig. 11 shows the variation of the peak frequency, Wp , with the peak 
attenuation, ilf p , for the various damping ratios f. For f — 0.1, the 
peak frequency is approximately the undamped frequency of the system 
and is independent of the peak attenuation. As the damping is increased, 
the peak frequency decreases, and it is lower for the lower peak attenua- 
tions. In general, the lower the value of Wp , the lower will be the band- 
width of the closed-loop system. 

The relation between the velocity gain constant wo and the peak 
magnitude is shown in Fig. 12 as a fimction of the damping ratio, K- 
As the gain constant is increased, the peak magnitude increases; in 
most cases (for constant peak magnitude), coo is less for the lightly 
damped systems. This graph shows that, for f = 0.1, a change in wo of 
approximately 3 db is sufficient to cause Mp to increase 8 db, while. 
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FREQUENCY RATIO 

Fig. S — Opcii-loop phase angle of G(s). 
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Fig. 9 — Gain margin va. closed-loop peak magnitude as a function of 
damping ratio. 
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Fig. 10 — Phase margin va, closed-loop peak magnitude as a function of damp- 
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Fig. 11 — Closed-loop peak frequency vs. peak magnitude as a function of 
damping ratio. 
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CLOSED LOOP PEAK MAGNITUDE IN DECIBELS 

Fig. 12 — (iiiiu ratio vs. cloHcd-loop peak magnitude as a functiou of damping 

ratio. 

for f = 1.00, wi, must change 9 db for the same increase m peak magni- 
tude. Since, in a practical system, fluctuations in the value of wu are to 
be expected, operation of Ihe lightly damped system would be more 
erratic tliau that of a .sy.stom having adequate damping. 

Equation (70), when solved for the transient response {di a step 
function) yield-s solutions as given in Fig. Ki. It is interesting to note 
that the transient response f(jr a servomechanisni having a damping 
ratio of ^ ^ 0.1 and a closed-loop peak magnitude Mp ^ 'S db shows 
little o\'erslioat. Examination of the frequency response shows that 
this is a result of a large attenuation in the frequency region below 
resonance. The superimjjosed o.scillation is caused by the gain in the 
resonant frequency region. 

The results of the transient solutions for t = 0-5 are summarized in 
graphical form in Fig. 14, in which the delay time, I'd, rise time, 1% , 
peak time, Tp , and per cent overshoot are given as functions of the 
gain, oju/c^n ■ (The definitions of the various time values are given in 
Fig. 15.) The response times decrease rapidly with increasing gain for 
the lower gain values and, as the gain increases, become relatively 
insensiti\e to gain variations. The per cent overshoot is a linear function 
of gain for values of gain above the limiting case in which there is no 
overshoot (oju/cij,, = O.'Mi^). 
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Fig. 13 — Variation of step response of linear sj'stem with damping ratio for 
constant closed-loop peak attenuation Mp = 3 db. 

The relationships developed in this section can be used to estimate 
design parameterB. For example, from Fig. 9, assume that a gain margin 
of 6 db is desired for a servomechanism having a damping ratio f = 0.5. 
This fixes the peak magnitude il/p to be 4.15 db and, from Fig. 10, the 
phase margin is found to be 50 degrees. From Fig. 11, the peak frequency 
is Wp = 0.8 w„ and, from Fig. 12, the gain ratio ojo/cj,. is —6 db; i.e., 
wo = 0.5 w„ . Fig. 14 then predicts an overshoot of 24 per cent; a delay 
time, 7'd = 2.3/a)„ ; a rise time, Tr = 2.1/w„ ; and a peak time, T-p = 
5.a/a;„ . 



VI. ANALOG SOLUTIONS OF THE NONLINEAR EQUATIONS 

The equations representing the behavior of a valve-controlled servo- 
mechanism were derived in Sections II and III, and the approximate 
linear theory was discussed in Sections IV and V. In this section, repre- 
sentative analog computer solutions of the nonlinear equations (which 
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Fig. 14 — Trausient response data for a damping riitio of 0.5. 

are suniniarized in Talile III) are given for particular values of the 
parameters, and the results are correlated with the linear solutions. 

The valve-controlled servomechani-sm is assumed to have the numeri- 
cal constants listed in Table I\'. During the course of the discussion, 
the effects of changes in these parameters will be considered, but, unless 
otherwise stated, the values will be assumed to be as given in the table. 
In this manner, a reference system is obtained and the discussion of 
the effects of parameter variations is facilitated by comparison with 
the reference l^ehavior. 

The first eight constants listed in the Table IV are considered to lie 
the independent variables while the remaining five are dependent. The 
compliance coefficient, /v, , is given by (IG) as the ratio of W to &. 
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The effective leakage coefficient, La, is given by (50) as 

Kbfo 



La — Lm + 



V2p, 



and the undamped natural frequency is given by (44) as 



(71) 



(72) 



Table III — Equations Solved by Analog Computation 
Qi = ■~L„(pi - Pi) - V^Bn + KrP2 

iPi - P'.)V„ = Jeo 

€ = 8i - do 
pi ^ and p« ^ 

Q, = Qm - Qei 

Q2 = Q02 - Qb«. 

A't I e + en I Vl P* — Pi | Sgn (p, — pi) e > — «o 
-0 t < -to 

"0 « > ta 

\_Kb I e — fo I Vl Ps- p2 I sgn (p. - pa) t < sq 
'0 « > (0 

]_J^b I e — to I Vpi e < eo 

A't I 6 + to i V^ e > — «0 

_0 e < -«o 



Qo2 = 
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Table I\' ^ — ^ Reference ^ 


ALUES OF PaR^VMETERS 


Definition 


Symbol 


Value 


Units 


total inertia 


J 


2.73 X 10-3 


lb -in-sec* 


motor (lisphieement 


F„ 


0.0151 


in^ 


trapped volume 
bulk modukis 


Vt 


0.125 


in" 


& 


2.22 X 10^ 


Ib/in^ 


motor leakage coefficient 


L,„ 


0.039 X 10-3 


in^/lb-sec 


supply pressure 


p.- 


1000 


Ib/in^ 


gain constant 


Kt 


0.0912 


in-Vlb"*-sce-rad 


open-center constant 


«o 


0.0561 


radians 






3.21 


degrees 


compliance coefficient 


K, 


0.0563 X 10-5 


inVlb 


effective leakage coefficient 


U 


0.1533 X 10-= 


in*/lb-sec 


resonant frequency 


Lin 


544.7 


rad/sec 






86.6 


cps 


damping ratio 


r» 


0.5 




gain ratio 




0.496 





The dimensionless damping ratio U [fi"om (51)] is 



fa y 



JP 



y 27 



and the velocity gain constant is obtained from (46) : 



Wa 






73) 



(74) 



Fig. l(i shows the theoretical frequency response of the reference 

servomechanism as a function of the hiput amphtude. The linear 
prediction (based on small amplitudes and pressure differentials) is 
included for comparison. The amplitude sensitivity in the small signal 
region, as represented for example by a curve of 1° amphtude, is the 
result of the nonhnear flow characteristics of the valve. The operation 
is within the region in which | e | < eo ( A region) and pressure saturation 
has not occurred, f 

The response for an input of 2° shows more deviation from the linear 
response, primarily because of pressure saturation; operation is still 
within the A region. For greater hiput ampHtudes, the response falls off 



t A —12 db per octave slope that passes through zero db at the frequency 



[/li. 



i)/ 



Je.- 



divides the graph into regions that represent the saturating and nonsaturating 
conditions. 
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Fig. 16 — Theoretical frequency response of a hj'draulic servomeclianism as a 
function of input amplitude. 



rapidly and operation is primarily in the B region; pressure saturation 
effects become more pronounced. 

The closed-loop transient response of the reference system as a func- 
tkm of step magnitude is given in Fig. 17, together with the linear 
solution. For small amplitudes (less than 3°), the linear and nonlinear 
solutions are es.sentially identical; as the amplitude is increased, the 
discrepancy becomes largo. The per cent overshoot decreases with 
amplitude. 

It is evident that the transient response of the servomechanism is 
not as sensitive to amplitude as is the frequency response; for example, 
comparison of the 3° curves in Figs. 16 and 17 shows that the transient 
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Fig. 17 — Theoretical transient response of hydraulic Rervomechanism as a 
function of amplitude. 

is much closer to the linear transient solution than the frequency re- 
sponse is to the linear frequency response. This difference is attributable 
mainly to pressure saturation. Fig. 18 shows the transient response of 
the system for the 3° step, together with the pressures pi and p2 . Since 
the supply pressure is 1000 psi and the valve is symmetrical, both 
pressures are initially 500 psi in the quie.scent condition. For the case 
shown, the applied .step was in the positive direction, so that ^o was 
also positive; consequently, pi has an initial positive slope while pi has 
a negative slope. From to 2 milliseconds, the oil is compressed and 
very little shaft rotation occurs; maximum acceleration (pi — pi) 
occurs at 2 milliseconds, at which time the shaft has acquired an ap- 
preciable velocity. From 2 to 5 milliseconds, the acceleration decreases 
from maximum to zero, while the velocity continues to increase to its 
maximum value. From 5 to 6.8 milliseconds, the acceleration becomes 
negative, since pi is now greater than pi . The velocity is still positive 
for this period, and the error signal e decreases from a positive value 
to zero. The increase in p2 is due to the compression of the oil in hne 2 
by the moving inertia, as the orifice area opening to the sump is gradu- 
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Fig. 18 — Output and pressure variation with time for a 3° step. 

ally reduced and the pump pressure is applied to this side of the motor. 
For the period from 6.8 to 9.5 milliseconds, the error e is negative and 
the velocity decreases to zero at the peak of the curve. 

The nature of the response after 9.5 milliseconds is very similar to 
that previously described, since the output starts from rest with an 
initial error. The major difference, aside from the fact that the error is 
now negative, is that the pressures have appreciable values at 9.5 
milliseconds, while at milliseconds the pressures started from the 
quiescent state. 

It should be noticed that the pressures pi and p2 did not limit for 
the 3° step and that, although the pressure differentials were appreciable, 
the linear theoiy still provided a good approximation to the output 
motion, as indicated in Fig. 17. The value of cq as given in Table IV is 
3.2°, so that operation was entirely within the A region. 

Fig. 19 shows the transient response and pressures for the 20° step. 
In this case, the pressures just limit at the start, and the pressure differ- 
entials are large. In Fig. 20, for the 30° step, the initial pressure satura- 
tion is more pronounced and, in addition, pi reaches zero and p2 obtains 
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Fig. 19 — Output ami pressure variation with time for a 20° step. 
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Fig. 20 — Output iiud pressure variation with time for a 30° step. 
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Fig. 21 — Output iiiid pressure variation with time for a 50° step. 



a value above that of the supply pressure (p„ = 1000 psi). Examination 
of the curves shows that the maximum value of p2 of 1080 psi occiu's 
for e positive, the velocity positive and the acceleration negative. This 
indicates that the inertia of the motor and load combination is forcing 
oil to flow out of the exhaust orifice and that, as a result, the pressure 
Pi achieves a very high value. At the same time, oil is forced through 
the orifice connecting the pressure supply to line 1, but the velocity is 
so great that the rate of flow is not sufficient to maintain a pressure in 
this line. Thus, during the period from 14 to 17 milliseconds cavitation 
conditions exist on this side of the motor. 

In Fig. 21, for the 50° step, the situation is similar to that of Fig. 20, 
except that the saturation and cavitation periods are of longer duration. 
Here, the peak pressure is very nearly 1400 psi. 

From the preceding discussion, the desirability of including relief 
valves in each line is apparent. Dangerously high pressures can be 
generated, especially if a supply pressure of 3000 psi is used. The inclu- 
sion of relief valv<'s in effect limits the line pressure, so that the oil 
cannot be "trapped" liy the inertia. This, however, has the disadvantage 
of increasing the overshoot for large amplitudes and does not decrease 
the rise or delay times. In addition, the cavitation period is prolonged. 

The effect of the variation in gain, Kb , is shown in Figs. 22 and 23. 
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Fig. 22 — Theoretical traneient response as a function of amplitude for Ki = 0.12. 

In Fig. 22, JU has a value of 0.12, f^o that the parameters listed on the 
figure change from those of the reference ease; all other parameters 
remain constant as given in Table IV. Since the open-center constant 
eo is now only 2.44°, it is to 1)0 expected that the system exhibit more 
amphtude sensitivity. Comparison with Fig. 17 shows that this is the 
case; in the small-signal region, the curves for the higher-gain system 
differ somewhat more from the linear solution than do those in Fig. 17. 
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Fig. 23 — Theoretical transient response as a function of amplitude for Kb = 0.184. 



154G THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1959 



<t|a>l.4 



1.2 
1.0 

o.e 

0.6 
0.4 
0,2 





































r 


\ 
























— 


~1 


L-\ 


-\ 




f 


\ 


— 


/' 








/ 


^ 




\ 


/ 


\, 


— 


7 


— 


._j 


/ 


\ 


J 




\. 


/ 


s 


y. 


— 


:^ 


\^> 
















/ 






























/ 


















eo=\.\T' 






J 


' 






















/ 





























10 15 20 Z5 30 35 40 45 50 55 60 65 70 75 
TIME IN MILLISECONDS 

Fig. 24 — Transient response for ii 10° step with Kb = 0.25. 

As the amplitude increases, this effect is lost, and the major difference 
is found in the per cent overshoot. 

Ill Fig. 23, the gain constant, Kh , is 0.184. For this vahie the dimcn- 
sionless gain ratio is unity and the linear theory predicts zero gahi 
margin. There is a distinct difference in the nature of the response as the 
amplitude increases. For the \° step, the operation is entirely within 
the A region («o = 1.59°) and the oscillations are almost continuous. 
As the amplitude increases, the per cent overshoot decreases, since 
pressure saturation occurs. 

This effect on the stability is more pronounced at the higher gain 
vahies. In Fig. 24, a 10° stop is .shown for the case in which Kb = 0.25 
and the gain ratio is l.'M"). From the linear theory, this system should 
be unstable. Examination of the response shows that, in the A region 
in which c is less than eo = 1.17°, the system is unstable, but that it is 
stable in the B region in which « is greater than to . The response, there- 
fore, oscillates indefinitely, but only with the amplitude of eo ■ It is 
evident, therefore, that, if frictional forces are sufficient to overcome 
the small oscillations of ±to amplitude, or if these oscillations are not 
detrimental to the performance in the particular apphcation, the allow- 
able gain is much greater than that predicted by the linear theory. For 
a given valve, /Ci,eo is a constant, so that, as the gain K), is increased, 
the magnitude co of the sustained oscillations decrease. The incremental 
damping and gain constants for the B region are given by (57) and (59) ; 
these equations show that the operation in the B region is inherently 
more stable than that in the A region, since the incremental damping 
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Fig. 25 — TraiiHieiit response as a function of supply pressure for a 20° step. 

is greater than in the A region and the incremental gain is less than the 
corresponding A region gain. Thus, in the absence of friction or viscous 
damping, the sustained oscillations are most pronounced in the A region, 

The effect of the supply pressure on the transient response is shown 
in Fig. 25 for a 20° step. The linear theory provides the best approxima- 
tion for the case in which the supply pressure is greatest. This is the 
result of two factors: (a) the system with the higher pressure is less 
susceptible to pressure saturation and (b) the value of co increases with 
increasing supply pressure, so that operation is more completely in the 
A region. The solutions show that, whereas the system having a supply 
pressure of 1000 psi encounters saturation for a 20° step, the 3000-psi 
servo is not pressure-limited until subjected to a 60° step. The advantage 
of operating at higher pressures is thus primarily a question of pressure 
saturation. 

In all the previous solutions, the effective leakage coefficient. La , 
has been maintained constant at 0.1533 X 10"^ inches^ per pound- 
second, so that the damping ratio, f„ , was 0.5. From (71), it is seen 
that the effective leakage coefficient is the sum of the motor leakage 
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Fig. 26 — Theoretical transient response as a function of Kitts for a 1° step. 

coefficient and a term proportional to Kh^a . It is interesting to observe 
the response of the system as affected by the nature of the damping. 
This is given in Fig. 2G for a 1° step. The three upper curves have a 
damping ratio of fa = 0.5; the difference in the curves results from the 
method used in obtaining the damping. For cases in which the open- 
center valve provides appreciable damping, the response does not differ 
from the linear prediction. This is the result of operation in the A region, 
where e is less than en . As the valve damping decreases and the motor 
leakage is increased, the response is slower and falls below^ the linear 
curve; this occurs for the case in which Kb€D is 0.0005. Operation is 
partly in the B region, since to is equal to 0.31". 

When the damping is contributed entirely by the motor leakage, 
Kbeo = and the valve is of the critical-center type. In this case. Fig. 26 
shows that the response does not overshoot and that an increase in 
gain would be desirable. Computer results show that, for Kb = 0.20, 
the critical-center valve gives a transient response having about 25 per 
cent overshoot and adequate stabihty. However, it should be empha- 
sized that the damping ratio for this case was 0.5 and that the damping 



VALVE-CONTROLLED HYDRATTLIC SERVOMECH.A.NISMS 1549 

was obtained by altering the motor leakage. Any attempt to reduce the 
valve-quiescent flow without compensating the system in some manner 
to provide additional damping results in an miderdamped servomecha- 

nism. 

The final curve, for which /u^o = -0.001 in Fig. 26, represents the 
response to be expected with a closed-center valve. This valve has a 
dead zone of 0.03° about which the output will wander, and is shown 
in a particularly poor case, since the step is only 1° and the gain is small. 
The only damping in this system is contributed by motor leakage, so 
that the system is underdamped. 

V!I. rONCLUSIOXS 

This study shows that the linear approximation to the nonlinear 
representation of valve-controlled hydraulic servomechanisms can be 
applied only with the sacrifice of considerable accuracy. However, since 
the linear theory is readily applied, it can be used in obtaining estimates 
for preliminary designs if the deviations from the nonlinear solutions 
are understood. 

REFERENOKS 

1 Schiirfer, J, W., An Electrically Operated Hydraulic Control Valve, B.R.T.J., 
3G, May lOf)?, p. 711. " . , ,. ^, , , 

2. Blarkbuni, J. F. and Lee, S. Y., Contributions to Hydraulic Coulrol — I. 
Steady-Slate Axial Forces on Control-Valve Pistons, Trans. A.S.M.E., 74, 
AuKUSt 1952, p. 1005, , . „ 

3 Blarldnirn, J. F., Contributions to Hydraulic Control — 3. Pressure-Flow 

Kelalionship.'s for 4-Wav Valves, Trans. A.H.M.E., 75, August 1953, p. 1163. 

4 Shearer J L., JJvnamic Characteristics of Valve-Controlled Hydraulic Servo- 

motors, Trans." A.S.M.E., 76, August 1954, p. 895. 



